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' Abstract 

We consider a variant of two-point Euclidean shortest path query problem: given a 
\Q • polygonal domain, build a data structure for two-point shortest path query, provided that 

04 | query points always lie on the boundary of the domain. As a main result, we show that 

a logarithmic-time query for shortest paths between boundary points can be performed 
\ using 0(n 5 ) preprocessing time and 0(n 5 ) space where n is the number of corners of the 

r \ • polygonal domain and the O- notation suppresses the poly logarithmic factor. This is realized 

by observing a connection between Davenport-Schinzel sequences and our problem in the 
parameterized space. We also provide a tradeoff between space and query time; a sublinear 
time query is possible using 0(n 3+e ) space. Our approach also extends to the case where 
query points should lie on a given set of line segments. 

> : 

■ 1 Introduction 

O : 

' A polygonal domains V with n corners and h holes is a polygonal region of genus h whose 

boundary consists of n line segments. The holes and the outer boundary of V are regarded as 
obstacles. Then, the geodesic distance between any two points p, q in a given polygonal domain 
^ \ V is defined to be the length of a shortest obstacle-avoiding path between p and q. 

The Euclidean shortest path problem in a polygonal domain has drawn much attention 
in the history of computational geometry [15]. In the two-point shortest path query problem, 
we preprocess V so that we can determine a shortest path (or its length) quickly for a given 
pair of query points p,q £ V . While we can compute a shortest path in 0(n log n) time from 
scratch [13], known structures for logarithmic time query require significantly large storage [6]. 
Chiang and Mitchell [6] developed several data structures that can answer a two-point query 
quickly with tradeoffs between storage usage and query time. Most notably, O(logn) query 
time can be achieved by using 0(n n ) space and preprocessing time; sublinear query time by 
0(n 5+e ) space and preprocessing time. More recently, Guo et al. [10] have shown that a data 
structure of size 0(n 2 ) can be constructed in 0(n 2 logn) time to answer the query in 0(h log n) 
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Table 1: Summary of new and known results on exact two-point shortest path queries, where 
e > is arbitrary and < 5 < 1 is a parameter, [new] denotes our results. 



Query domain 


Preprocessing time 


Space 


Query time 


Ref. 


V 


0(re n ) 


0(re n ) 


O(logre) 


[6] 


V 


0(n 10 logn) 


0(n 10 log re) 


0(log 2 re) 


[6] 


V 


O(n 5+105+e ) 


0(n 5+w5+e ) 


0(re 1_<5 log re) 


[6] 


V 


<3(n 5 ) 


0(re 5 ) 


0(logre + h) 


[6] 


V 


0(re + /i 5 ) 


0{n + h 5 ) 


0(h log re) 


[6] 


V 


0(n 2 log n) 


0{n 2 ) 


0(h log re) 


[10] 


&P 


0(re 4 A 65 (re)logn) 


0(n 4 A 66 (re)) 


O(logre) 


[new] 


&P 


0(n 3+s X 65 {n 5 )logn) 


0(n 3+5 X m (n s )) 


0(re 1-5 log re) 


[new] 


m segments 


0(m 2 n 3+5 X^ (n 5 ) log re) 


0(m 2 n 3+5 X m (n s )) 


0(re 1_<5 log(m + re)) 


[new] 



time, where h is the number of holes. Their results are summarized in Table [TJ For more results 
on shortest paths in a polygonal domain, we refer to survey articles by Mitchell [14, 15]. 

In this paper, we focus on a variant of the problem, in which possible query points are 
restricted to a subset of V\ the boundary dV of the domain V or a set of line segments within 
V . In many applications, possible pairs of source and destination do not span the whole domain 
V but a specified subset of V . For example, in an urban planning problem, the obstacles 
correspond to the residential areas and the free space corresponds to the walking corridors. 
Then, the query points are restricted to the spots where people depart and arrive, which are on 
the boundary of obstacles. 

Therefore, our goal is to design a data structure using much less resources than structures 
of Chiang and Mitchell [6] when the query domain is restricted to the boundary of a given 
polygonal domain V or to a set of segments in V . To our best knowledge, no prior work 
seems to investigate this variation. As a main result, in Section [3J we present a data structure 
of size 0(re 4 A 66 (n)) that can be constructed in 0{n X^{n) log re) time and can answer a dV- 
restricted two-point shortest path query in O(logre) time. Here, X m (n) stands for the maximum 
length of a Davenport-Schinzel sequence of order m on re symbols [19]. It is good to note 
that A m (re) = 0(n log* re) for any constant m as a convenient intuition, while tighter bounds 
are known [16, 19]. We also provide a tradeoff between space and query time in Section [TJ 
In particular, we show that one can achieve sublinear query time using 0(re 3+<E ) space and 
preprocessing time. New results in this paper are also summarized in Table [TJ 

Our data structure is a subdivision of two-dimensional domain parameterized in a certain 
way. The domain is divided into a number of grid cells in which a set of constrained shortest 
paths between query points have the same structure. Each grid cell is divided according to the 
projection of the lower envelope of functions stemming from the constrained shortest paths. 
With careful investigation into this lower envelope, we show the claimed upper bounds. 

Also, our approach readily extends to the variant where query points are restricted to lie on 
a given segment or a given set of segments in V . We discuss this extension in Section 

1.1 Related Work 

In the case where V is a simple polygon {h = 0), the two-point shortest path query can be 
answered in O(logre) time after 0{n) preprocessing time [9]. More references and results on 
shortest paths in simple polygons can be found in a survey article by O'Rourke and Suri [17] 

Before Chaing and Mitchell [6], fast two-point shortest path queries in polygonal domains 
were considered as a challenge. Due to this difficulty, many researchers have focused on the 
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approximate two-point shortest path query problem. Chen [5] achieved an 0(n log n)-sized 
structure for (6 + e)-approximate shortest path queries in O(logn) time, and also pointed out 
that a method of Clarkson [7] can be applied to answer (l+e)-approximate shortest path queries 
in O(logn) time using 0(n 2 ) space and 0(n 2 log n) preprocessing time. Later, Arikati et al. [4] 
have improved the above results based on planar spanners. 

The problem on polyhedral surfaces also have been considered. Agarwal et al. [1] presented 
a data structure of size 0(n 6 m 1+ ) that answers a two-point shortest path query on a given 
convex polytope in O {{^Jn/m l I A ) log n) time after 0(n 6 m 1+s ) preprocessing time for any fixed 

1 < m < n 2 and any 5 > 0. They also considered the problem where the query points are 
restricted to lie on the edges of the polytope, reducing the bounds by a factor of n from the 
general case. Recently, Cook IV and Wenk [8] presented an improved method using kinetic 
Voronoi diagrams. 

2 Preliminaries 

We are given as input a polygonal domain V with h holes and n corners. More precisely, V 
consists of an outer simple polygon in the plane M 2 and a set of h (> 0) disjoint open simple 
polygons inside P. As a set, V is the region contained in its outer polygon excluding the holes, 
also called the free space. The complement of V in the plane is regarded as obstacles so that any 
feasible path does not cross the boundary &P and lies inside V . It is well known from earlier 
works that there exists a shortest (obstacle- avoiding) path between any two points p,q 6? [14]. 

Let V be the set of all corners of V . Then any shortest path from p £ P to g 6 P is a 
simple polygonal path and can be represented by a sequence of line segments connecting points 
in V U {p, q} [14]. The length of a shortest path is the sum of the Euclidean lengths of its 
segments. The geodesic distance, denoted by d(p,q), is the length of a shortest path between p 
and q. Also, we denote by \pq\ the Euclidean length of segment pq. 

A two-point shortest path query is given as a pair of points (p, q) with p,q E V and asks 
to find a shortest path between p and q. In this paper, we deal with a restriction where the 
queried points p and q lie on the boundary dV. 

A shortest path tree SPT(p) for a given source point p E V is a spanning tree on the corners 
V plus the source p such that the unique path to any corner v G V from the source p in SPT{p) 
is a shortest path between p and v. The combinatorial complexity of SPT{p) for any p G V is 
at most linear in n. A shortest path map SPM(p) for the source p is a decomposition of the 
free space V into cells in which any point x has a shortest path to p through the same sequence 
of corners in V. Once SPT(p) is obtained, SPM(p) can be computed as an additively weighted 
Voronoi diagram of V U {p} with weight assigned by the geodesic distance to p [14]; thus, the 
combinatorial complexity of SPM(p) is linear. A cell of SPM(p) containing a point q G V has 
the common last corner v G V along the shortest path from p to q; we call such a corner v 
the root of the cell or of q with respect to p. An 0(n log n) time algorithm, using O(nlogn) 
working space, to construct SPT(p) and SPM{jp) is presented by Hershberger and Suri [13]. 

An SPT- equivalence decomposition A SPT of V is the subdivision of V into cells in which every 
point has topologically equivalent shortest path tree. An A SPT can be obtained by overlaying 
n shortest path maps SPM(v) for every corner v G V [6]. Hence, the complexity of A SPT is 
0(n 4 ). Note that A SPT n dV consists of at most 0(n 2 ) points; they are intersection points 
between any edge of SPM(v) for any v G V and the boundary dV . We call those intersection 
points, including the corners V, the breakpoints. The breakpoints induce 0(n 2 ) intervals along 
dV. We shall say that a breakpoint is induced by SPM(v) if it is an intersection of an edge of 
SPM(v) and dV. 

Given a set V of algebraic surfaces and surface patches in R d , the lower envelope £(T) of 
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r is the set of pointwise minima of all given surfaces or patches in the d-th. coordinate. The 
minimization diagram M(T) of T is a decomposition of IR^" 1 into faces, which are maximally 
connected region over which C(F) is attained by the same set of functions. In particular, when 
d = 3, the minimization diagram M{T) is simply a projection of the lower envelope onto the 
xy-plane. Analogously, we can define the upper envelope and the maximization diagram. 

As we intensively exploit known algorithms on algebraic surfaces or surface patches and 
their lower envelopes, we assume a model of computation in which several primitive operations 
dealing with a constant number of given surfaces can be performed in constant time: testing if 
a point lies above, on or below a given surface, computing the intersection of two or three given 
surfaces, projecting down a given surface, and so on. Such a model of computation has been 
adopted in many research papers; see [2,3,18,19]. 

3 Structures for Logarithmic Time Query 

In this section, we present a data structure that answers a two-point query restricted on &P in 
O(logn) time. To ease discussion, we parameterize the boundary dV. Since dV is a union of 
h+1 closed curves, it can be done by parameterizing each curve by arc length and merging them 
into one. Thus, we have a bijection p: [0, \ &P\) — ► dV that maps a one-dimensional interval into 
dV, where \dV\ denotes the total lengths of the h + 1 closed curves forming &P . Conversely, 
the inverse of p maps each interval along dV to an interval of [0, 

A shortest path between two points p, q G V is either the segment pq or a polygonal chain 
through corners in V. Thus, unless d(p, q) = \pq\, the geodesic distance is taken as the minimum 
of the following functions f UjV : [0, \&P\) x [0, \&P\) — > M over all u,v G V, which are defined as 
follows: 



where VP(x), for any point x £ V, denotes the visibility profile of x, defined as the set of all 
points y G V that are visible from x; that is, xy lies inside V . The symbol oo can be replaced 
by an upper bound of max Si< d(p(s),p(t)); for example, the total length \dV\ of the boundary 
of the polygonal domain V. 

Since the case where p(s) is visible from p(t), so the shortest path between them is just the 
segment p(s)p(t), can be checked in O(logn) time using 0(n 2 logn) space [6], we assume from 
now on that p(s) ^ VP(p{t)). Hence, our task is to efficiently compute the lower envelope of 
the 0(n 2 ) functions f UjV on a 2-dimensional domain V := [0, \ &P\) x [0, \ &P\). 

3.1 Simple lifting to 3-dimension 

Using known results on the lower envelope of the algebraic surfaces in 3-dimension, we can show 
that a data structure of size 0(n 6+e ) for O(logn) query can be built in 0(n 6+e ) time as follows. 

Fix a pair of intervals I s and It induced by the breakpoints. Since I s belongs to a cell of an 
SPT-equivalence decomposition, VP(p(s)) is independent of choices over all s G I s . Therefore, 
the set V s := V H VP(p(s)) of corners visible from p(s) is also independent of the choice of 
s G I s and further, for a fixed u G V s , there exists a unique v G V that minimizes f u ,v(s,t) for 
any (s,t) G I s x It over all v G V [6]. This implies that for each such subdomain I s x It C T> 
we extract at most n functions, possibly appearing at the lower envelope. Moreover, in I s x I t , 
such a function is represented explicitly; for u G V s and v G Vt, 




\p(s)u\ + d(u,v) + \vp{t)\ if u G VP(p(s)) and v G VP{p{t)), 
co otherwise, 




(s, t) = ^/{x{s)-x u ) 2 + (y(s)-y u ) 2 + d(u, v) + ^{x{t)-x v ) 2 + (y(t)-y v ) 2 , 
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where x(s) and y(s) are the x- and the y-coordinates of p(s), and x u and y u are the x- and 
the y-coordinates of a point u G M 2 . Note that x(s) and y(s) are linear functions in s by our 
parametrization. 

For each u G V, there exists a corner v £ Vt that minimizes in I s x I t over all v' G Vj. 
Therefore, the function ^ u := min^y on I s x is well-defined. Observe that the graph of 
g u is an algebraic surface with degree at most 4 in 3-dimensional space. Applying any efficient 
algorithm that computes the lower envelope of algebraic surfaces in M 3 , we can compute the 
lower envelope of the functions g u in 0(n 2+e ) time [18]. Repeating this for every such subdomain 
I s x It yields 0(n 6+e ) space and preprocessing time. 

Since we would like to provide a point location structure in domain T>, we need to find 
the minimization diagram Ai of the computed lower envelope. Fortunately, our domain is 2- 
dimensional, so we can easily project it down on V and build a point location structure with 
an additional logarithmic factor. 

In another way around, one could try to deal with surface patches on the whole domain T>. 
Consider a fixed corner u G V and its shortest path map SPM(u). The number of breakpoints 
induced by SPM(u) is at most O(n), including the corners V themselves. This implies at 
most an 0(n 2 ) number of combinatorially different paths between any two boundary points 
p(s) and p(t) via u. That is, for a pair of intervals I s and It, we have a unique path via u 
and its length is represented by a partial function of (s, t) defined on a rectangular subdomain 
I s x I t C V. Hence, we have 0(n 2 ) such partial functions for each u G V, and thus 0(n 3 ) 
in total. Each of them defines an algebraic surface patch of constant degree on a rectangular 
subdomain. Consequently, we can apply the same algorithm as above to compute the lower 
envelope of those patches in 0((n 3 ) 2+e ) = 0(n 6+e ) time and space. 

3.2 0(n 5+e )-space structure 

Now, we present a way of proper grouping of subdomains to reduce the time/space bound by 
a factor of n. We call a subdomain I s x I t C V, where both I s and It are intervals induced by 
breakpoints, a grid cell. Thus, V consists of 0(n 4 ) grid cells. We will decompose V into 0(n 3 ) 
blocks of 0(n) grid cells. 

Consider a pair of boundary edges S,T C dV and let bs and br be the number of breakpoints 
on S and on T, respectively. Let p(sq), . . . ,p(sb s ) and p(to), . . . ,p(tb T ) be the breakpoints on S 
and on T, respectively, in order so < si < • • • < Sb s and to < t± < ■ ■ ■ < tf, T - Take br grid cells 
with s G [so) 5 !.) and let C := [sq,si) x [to,tb T ) C [so,S{, s ) x [to,tf, T ) be their union. We redefine 
the functions f U:V on domain C. As discussed above, for any s G [so,si), we have a common 
subset V s of corners visible from p(s). 

For any u G V s , let g u (s,t) := min^gy fu,v(s, t) be a function defined on C and b^ be the 
number of breakpoints on T induced by SPM(u). The following is our key observation. 

Lemma 1 The graph of g u {s, t) on C consists of at most 6^ + 1 algebraic surface patches with 
constant maximum degree. 

Proof. If g u (s,t) = f u ,v(s,t) for any (s,t) G C and some v G V, then pit) lies in a cell of 
SPM{u) with root v; by the definition of g u , the involved path goes directly from p(s) to u and 
follows a shortest path from u to On the other hand, when we walk along T as t increases 
from t\ to tb T , we encounter b^ breakpoints induced by SPM(u); thus, 6^ + 1 cells of SPM(u). 
Hence, the lemma is shown. ED 

Moreover, the partial function corresponding to each patch of 7 U is defined on a rectangular 
subdomain [sq, si) x [t^, ij) for some 1 < i < j < br- This implies that the lower envelope of g u 
on C is represented by that of at most ^2 u (b^ + 1) = n + br surface patches. 
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Though this envelope can be computed in 0((n + 6t) 2+€ ) time, we do further decompose C 
into [~^r-] blocks of at most n grid cells. This can be simply done by cutting C at t = t{ n for 
each i = 1, . . . , [^-J • For each such block of grid cells, we have at most In surface patches and 
thus their lower envelope can be computed in 0(n 2+e ) time. Hence, we obtain the following 
consequence. 

Theorem 1 One can preprocess a given polygonal domain V in 0(n 5+e ) time into a data 
structure of size 0(n 5+e ) that answers the two-point shortest path query restricted to the 
boundary dV in 0(logn)-time, where e is an arbitrarily small positive number. 

Proof. Recall that *£ s b s = ^2 T b T = 0(n 2 ). For a pair of boundary edges S and T, we can 
compute the lower envelope of the functions f U)V in 0(bs\^~\n 2+e ). Summing this over every 
pair of boundary edges, we have 

£ 0{b s \ h -^W^) = Oin^) ■ J2(~ + 1) = 0(n^). 
z — ' n z — ' n 

S,T T 

A point location structure on the minimization diagram can be built with additional logarithmic 
factor, which is subdued by 0(n e ). ED 



3.3 Further improvement 

The algorithms we described so far compute the lower envelope of surface patches in 3-space 
in order to obtain the minimization diagram Ai of the functions f uv . In this subsection, we 
introduce a way to compute M rather directly on T>, based on more careful analysis. 

Basically, we make use of the same scheme of partitioning the domain T> into blocks of (at 
most) n grid cells as in Section [3j 2. Let C be such a block defined as [so, s\) x [to, tb T ) such that 
[so,si) is an interval induced by the breakpoints and [toi4 T ) ^ s a union of bx < n consecutive 
intervals in which we have 6y — 1 breakpoints t\, . . . , 4 T -l- 

By Lemma Q3 the functions g u = min„ G y f u>v restricted to C can be split into at most In 
partial functions hi with 1 < i < 2n defined on a subdomain Cj C C. Each /ij(s, t) is represented 
explicitly as hi(s,t) = \p(s)ui\ + d(iii,Vi) + |fi£>(t)| in C{, so that we have hi(s,t) = f Ui ,Vi(s,t) 
for any (s,t) £ Ci and some Ui,Vi G V. Note that it may happen that Ui = uj or V{ = vj for 
some i and j; in particular, if Ui = Uj, we have Ci C\Cj = 0. Also, as discussed in Section [3J2, Ci 
is represented as [so^i) x [tk,tk') for some < k < k' < bx- 

In this section, we take the partial functions hi into account, and thus the goal is to compute 
the minimization diagram M of surface patches defined by the hi. We start with an ordering 
on the set V s of corners visible from p{s) for any s £ [sq, si) based on the following observation. 

Lemma 2 The angular order of corners in V s seen at s is constant if s varies within [sq, si). 

Proof. If this is not true, we have such an s' G [so, si) that p(s') and two corners v , v' in are 
collinear. Since corners lie on the boundary of an obstacle, one of v and v' is not visible from 
p(s) locally near s'; that is, p(s') is a breakpoint, a contradiction. ED 

Without loss of generality, we assume that as s increases, p(s) moves along dV in direction 
that the obstacle lies to the right; that is, p(s) moves clockwise around each hole and counter- 
clockwise around the outer boundary of V . (See Figure [TJ) By Lemma [21 we order the corners 
in in counter-clockwise order at p(s) for any s £ [sq,si); let -< be a total order on V s such 
that u -< v! if and only if Zp(so)p(s)u < Zp(so)p(s)u' . 

From now on, we investigate the set 

B(ij) -.= {(s,t) ednCj | hi(s,t) = hj(s,t)}, 



6 



Figure 1: As s increases, p(s) moves in direction where the obstacle is to the right. 



which is a projection of the intersection of two surface patches defined by h{ and hj. One can 
easily check that B(i,j) is a subset of an algebraic curve of degree at most 8. 

Lemma 3 The set B(i,j) is t-monotone. That is, for fixed t, there is at most one s £ [sq, s\) 
such that (s,t) £ B(i,j). 

Proof. If (s,t) £ B(i,j), we have equation hi(s,t) = hj(s,t). From the equation, we get 
\p(s)ui\ — \p(s)uj\ = \vjp(t)\ + d(uj,Vj) — \vip(t)\ — d(ui,Vi). If we fix t as constant, \p(s)ui\ — 
\p(s)uj\ remains a constant even if s varies within [so,si). Thus, p(s) is an intersection point 
with line segment p(so)p(si) and a branch H of a hyperbola whose foci are Ui and Uj. If 
p(so)p(si) C\H consists of two points, then the line £ through Ui and Uj must cross p(so)p(si) at 
a point z; I is the transverse axis of H and both Ui and Uj lie in one side of the line supporting 
p(so)p(si). Such a crossing point z is a breakpoint by definition but we do not have any break- 
point within p(so)p(si) since [so, si) is an interval induced by the breakpoints, a contradiction. ED 

Lemma 4 The set B{i,j) is either an empty set or an open curve whose endpoints lie on 
the boundary of Ci n Cj. Moreover, B(i,j) is either a linear segment parallel to the t-axis or 
s-monotone. 

Proof. Let I s and It be intervals such that Cj n Cj = I s X It- Note that I s = [so, si) and 
h = [tk,tk>) for some < k < k' < bx- 

First, note that if Ui = Uj, Ci C)Cj =0 and thus B(i,j) = 0, so the lemma is true. Thus, we 
assume that Ui ^ Uj. Regarding Vi and Vj, there are two cases: Vi = Vj or Vi ^ vj. In the former 
case, we get |p(s)uj| — \p(s)uj\ = d(v,j,Vj) — d(ui,Vi) from equation hi(s,t) = hj(s,t). Observe 
that variable t is readily eliminated from the equation, and thus if there exists (s',t') G Cj D Cj 
with (s',t') £ B(i,j), we have (s',t) £ B(i,j) for every other t £ It- Hence, by Lemma[3j B(i,j) 
is empty or a straight line segment in Ci n Cj which is parallel to t-axis, and thus the lemma is 
shown. 

Now, we consider the latter case where Vi ^ Vj. Without loss of generality, we assume that 
Ui ~< Uj. Recall that if U{ -< Uj, then Zp(so)p(s)ui < Zp(so)p(s)uj for any s in the interior of I s . 
We denote 0j(s) := Zp(so)p(s)ui and Oj(s) := Zp(so)p(s)uj . On the other hand, we also have 
a similar relation for Vi and Vj. Let (f>i(t) := Z.p(to)p(t)vi and <pj(t) := Zp{to)p{t)vj. Observe 
that 4>i(t) and (f>j(t) are continuous functions of t, and if 4>i(t') = 4>j(t') at t = t', then p(t') is 
a breakpoint induced by SPM{ui) or SPM(uj). Since It contains no such breakpoint induced 
by SPM(ui) or SPM(uj) in its interior, either </>i(t) < 4>j(t) or <pi{t) > 4>j(t) for all t in the 
interior of It] that is, the sign of 4>j(t) — 4>i{t) is constant. 

Since for any s,s' £ I s with s' > s we have \p(s')p(s)\ = s' — s by our parametrization, we 



can represent \p(s)ui\ = y (s + a^) 2 + bj and = y (t + q) 2 + df, where Oj, bi, Ci and di 

are constants depending on u^, Vi, and parametrization p. More specifically, s + a% denotes a 
signed distance between p(s) and the perpendicular foot of Uj onto the line supporting p(I s ), 
and bi is the distance between Uj and the line supporting p(I s ). See Figured Thus, hi(s,t) can 
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represented as hi(s, t) = y(s + a«) 2 + b\ + J (t + q) 2 + dj + d(ui,Vi). 

Now, we differentiate the both sides of equation hi(s, t) = hj(s, t) by t to obtain the derivative 

ds . 
dt ■ 

s + di ds t + Ci s + dj ds t + Cj 

y/( 8 + ai y + q dt + ft + Ci y + d\ ~ ft + a 3 y + b) dt + + Cj ) 2 + d) ' 

Rearranging this, we obtain 

t + C{ t + Cj 
7^^^^=^^^^= H 

ds _ )J( t + a) 2 + j y/( t + c i) 2 + d ) __ -cos(^(t)) +cos(^(t)) 
dt ~ g + Qj f_+ Qj ~ cos(^(s)) - cos(^(s)) 

^( s + ai )2 + ^ ^/( S + a i ) 2 + ^ 2 

Since < ^i(s) < Oj(s) < ir, we have cos(0j(s)) — cos(^(s)) > 0. Also, as discussed above, 
— (j>i(t) has a constant sign when t varies within the interior of It- Thus, cos(<j>j(t)) — 
cos(4>i(t)) has a constant sign, and ^ also has a constant sign at any (s,t) E B(i,j). Further- 
more, t| is continuous and has no singularity in the interior of C,, D Cj. This, together with 
Lemma [3l proves the lemma. EO 



Now, we know that B(i,j) can be seen as the graph of a partial function {s = j(t)}. Also, 
Lemma H] implies that B(i,j) bisects CiCiCj into two connected regions R(i,j) and R(j,i), where 
R(i,j) := {(s,t) e dnCj I hi(s,t) < hj(s,t)} and R(j,i) := {(s,t) G CiPlCj | hi(s,t) > hj(s,t)}. 
Let be the set of points (s,t) where the minimum of hj(s,t) over all j is attained by 

h{(s, t). We then have M(i) =Ci\ (Jj R(j, i) for each i. 

For easy explanation, from now on, we regard the s-axis as the vertical axis in T> so that we 
can say a point lies above or below a curve in this sense. 

The idea of computing M.(i) is to use the lower and the upper envelopes of the bisecting 
curves B(i,j). In order to do so, we extend B(i,j) to cover the whole t-interval It = [tk,tk') 
in C{ n Cj by following operation: For each endpoint of B(i,j), if it does not lie on the vertical 
line {t = tk} or {t = attach a horizontal segment to reach the vertical line as shown 

in Figure [3^a). We denote the resulting curve by (3(i,j); if B(i,j) = 0, define f3(i,j) as the 
horizontal segment connecting two points (so, tk) and (sojifc') i n CidCj. Observe now that (3(i,j) 
bisects Ci n Cj into regions R(i,j) and R(j,i), which lie above and below f3(i,j), respectively. 

Let P(i,j) + C dnCj be the region above f3(i,j) and (3(i,j)~ be the region below (3(i,j). For 
a fixed i with 1 < j < 2n, we classify the (3(i,j) into two sets Lj and Ui such that (3(i,j) £ if 
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Figure 3: (a) How to extend B(i,j) (dashed line) to (3(i,j) by attaching horizontal segments 
(solid line); (b) £(Lj), (c) U(Ui), and (d) the region M.(i) between them. Here, the dotted 
boxes are grid cells whose union is Cj and the s-axis appears vertical. 



R(j,i) = or 0{i,j) e Ui if R(j,i) = P(i,j)-. Recall that Mii) = d \ U, R(j,i). Thus, 

we have 

M(i)=Ci \ ( (J /3+U [J p- \ ■ 

The boundary of U/3eLi @ + m is t ne lower envelope C{L{) of Lf, symmetrically, the boundary 
of [_}p < z Ui P~ in Ci is the upper envelope U(Ui) of Ui. Therefore, M.{i) = C{L,i)~ nU(Ui) + , the 
region below the lower envelope C{Lj) of Li and above the upper envelope U{Ui) of Ui, and it 
can be obtained by computing the overlay of two envelopes C{Li) and U(Ui). See Figure [3](b) 
(d). We exploit known results on the Davenport-Schinzel sequences to obtain the following 
lemma [12,19]. 

Lemma 5 The set A4(i) is of combinatorial complexity 0(\QQ(n)) and can be computed in 
0(\65(n) logn) time, where A m (n) is the maximum length of a Davenport-Schinzel sequence of 
order m on n symbols. 

Proof. consists of at most three arcs, at most one algebraic curve of degree 8 and at 

most two straight segments. Thus, we have at most 6n algebraic arcs of degree at most 8. 
Any two such arcs can intersect each other at most 64 times by Bezout's Theorem [11]. Thus, 
each of C{Li) and U{Ui) has complexity 0(\§§{n)) and can be computed in 0{\^{n) log n) 
time [12,19]. 

After sorting the vertices on these envelopes in i-increasing order, we can easily specify all 
intersections between C{L.{) and U(Ui) in the same bound. ED 



It should be noted here that the exact constant 66 is not relevant; it only matters that this 
is some constant. 

We can compute the minimization diagram M. by computing each in 0{n\^{n) logn) 

time. In the same time bound, we can build a point location structure on M. Finally, we 
conclude our main theorem. 

Theorem 2 One can preprocess a given polygonal domain V in 0{n A \^{n) logn) time into a 
data structure of size 0(n A \§§{n)) that answers the two-point shortest path query restricted on 
the boundary dV in O(logn) time. 
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4 Tradeoffs Between Space and Query Time 



In this section, we provide a space/query-time tradeoff. We use the technique of partitioning 
V, which has been introduced in Chiang and Mitchell [6]. 

Let 5 be a positive number with < S < 1. We partition the corner set V into m = re 1_<5 
subsets V\, . . . ,V m of near equal size 0(n 5 ). For each such subset V{ of corners, we run the 
algorithm described above with little modification: We build the shortest path maps SPM(u) 
only for u £ Vi and care about only 0(n 1+s ) breakpoints induced by such SPM(u). Thus, we 
consider only the paths from p(s) via u £ Vi and v £ V to p(t), and thus 0{n 1+s ) functions f u , v 
for u EVi and v £ V. 

Since we deal with less number of functions, the cost of preprocessing reduces from O(re) 
to 0(n s ) at several places. We take blocks of 0(n s ) grid cells contained in T> and the number 
of such blocks is 0(n 2+s ). For each such block, we spend 0(n s X^(n s ) log re) time to construct 
a point location structure for the minimization map Mi of the functions. Iterating all such 
blocks, we get running time 0(n 2+2S X^(n s ) log n) for a part Vi of V. Repeating this for all such 
subsets Vi yields 0(n 3+s XQ5(n s ) log re) construction time. 

Each query is processed by a series of m point locations on every Mi, taking 0(m log re) = 
0(n 1_<5 logre) time. 

Theorem 3 Let 5 be a fixed parameter with < 8 < 1. Using 0(n 3+s \Q^(n s ) log re) time 
and 0(n 3+s XeQ(n s )) space, one can compute a data structure for 0(n 1-5 log re) -time two-point 
shortest path queries restricted on the boundary dV . 

Remark that when 5 = 1, we obtain Theorem [21 and that 0(n 3+e ) time and space is enough 
for sublinear time query. Note that if 0(n) time is allowed for processing each query, 0(n 2 ) 
space and 0(n 2 log n) preprocessing time is sufficient. 

5 Extensions to Segments-Restricted Queries 

Our approach easily extends to the two-point queries in which queried points are restricted to 
be on a given segment or a given polygonal chain lying in the free space V. 

Let S s and St be two sets of m s and rrij line segments, respectively, within V . In this section, 
we restrict a query pair (p, q) of points to lie on S s and St each. We will refer to this type of 
two-point query as a (S s , St) -restricted two-point query. As we did above, we take two segments 
S G S s and T £ St and let bs and bx be the number of breakpoints — the intersection points 
with an edge of SPM(v) for some v £ V — on S and T, respectively. Also, parameterize S 
and T as above so that we have two bijections p : [0, \S\] — > S and q : [0, \T\] — > T. 

Any path from a point on S leaves to one of the two sides of S. Thus, the idea of handling 
such a segment within the free space V is to consider two cases separately. Here, we regard S 
and T as directed segments in direction of movement of p(s) and q(t) as s and t increases, and 
consider only one case where paths leave S to its left side and arrive at T from its left side. The 
other cases are analogous. 

Then, the situation is almost identical to that we considered in Section [3J.2. For a pair of 
segments S and T, we can construct a query structure in 0(bs \^-] n\Q^{n) log n) time. Unfortu- 
nately, bs and br can be as large as 0(n 2 ), yielding the same time bound for (dV, c^P )-restricted 
two-point queries in the worst case. Thus, in the worst case, we need additional factor of m s mt 
as follows. 

Theorem 4 Let S s and St be two sets of m s and mt (possibly crossing) line segments, respec- 
tively, withinV, andS be a fixed parameter with < S < 1. Then, using 0(m s mtn 3+s \^(n s ) log re) 
time and 0(m s mtn 3+s Xmi^)) space, one can compute a data structure for 0(n 1_<5 log(re+m s + 
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mt))-time (S s , St) -restricted two-point queries. 

Remark that in practice we expect that the number of breakpoints b$ and Y2t is 
not so large as f2(n 2 ) that the preprocessing and required storage would be much less than the 
worst case bounds. 

6 Concluding Remarks 

In this paper, we posed the variation of the two-points query problem in polygonal domains 
where the query points are restricted in a specified subset of the free space. And we obtained 
significantly better bounds for the boundary-restricted two-point queries than for the general 
queries. 

Despite of its importance, the two-point shortest path query problem for the polygonal 
domains is not well understood. There is a huge gap (0(n) to 0(n n )) about logarithmic query 
between the simple polygon case and the general case but still the reason why we need such 
a large storage is still unclear. On the other hand, restriction on the query domain provides 
another possibility of narrowing the gap with several new open problems: (1) What is the 
right upper bound on the complexity of the lower envelope defined by the functions f U)V on 
the parameterized query domain? And what about any lower bound construction? (2) If the 
query domain is a simple 2-dimensional shape, such as a triangle, then can one achieve a better 
performance than the general results by Chiang and Mitchell? 

We would carefully conjecture that our upper bound (9(n 5 ) for logarithmic query could be 
improved to 0(n 4 ). Indeed, we have 0(n 4 ) grid cells on the parameterized query domain and 
whenever we cross their boundaries, changes in the involved functions are usually bounded by 
a constant amount. Thus, if one could find a clever way of updating the functions and their 
lower envelope, it would be possible to achieve an improved bound. 
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